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ABSTRACT 

A  system  for  the  analysis  of  geomagnetic  and  gradient 
noise  is  described.   The  nominal  176  kHz  Larmor  frequency 
output  of  an  optically  pumped  cesium  magnetometer  is  mixed 
to  a  frequency  in  the  1.2  kHz  range  and  FM  tape  recorded. 
The  recorded  output  is  applied  to  a  periodmeter,  the  digital 
output  of  which  is  input  to  the  XDS-9300  computer  for 
analysis  and  display.   The  absolute  accuracy  of  the  period 
measurement  is  .75   usee  with  a  resolution  of  .25   ysec. 
The  RMS  noise  of  the  tape  recorder  was  determined  to  be 
.667  ysec  or  in  the  operating  frequency  range  less  than 
1  Hz  or  .28  gamma.   Externally  introduced  field  disturbances 
as  small  as  .5  gamma  were  recorded  and  displayed,  and  the 
response  time  of  the  magnetometer  determined  to  be  less 
than  2.5  ms  ec . 
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I.   INTRODUCTION 

With  a  continuing  interest  in  magnetic  anomaly  detec- 
tion the  Naval  Postgraduate  School  has  studied  geomagnetic 
noise,  first  with  the  proton  precession  magnetometer  and 
later  with  the  optically  pumped  rubidium  type.   In  1972 
the  school  obtained  an  optically  pumped  cesium  magnetometer 
with  accessories  to  convert  it  to  a  gradiometer.   In  order 
to  take  full  advantage  of  the  inherent  high  frequency 
response  of  this  magnetometer  a  data  collection  system 
was  developed  to  observe  geomagnetic  and  gradient  noise. 
Figure  1.1  is  a  block  diagram  of  the  system. 

A  brief  introduction  to  the  earth's  magnetic  field 
and  other  measuring  instruments  is  presented  first  followed 
by  a  description  of  the  system  developed. 
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Figure  1.1   System  Block  Diagram 
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II.   THE  MEASUREMENT  FIELD 

A.   MAIN  FIELD 

World  magnetic  surveys  both  at  ground  level  and  more 
recently  by  satellite  indicate  that  the  principal  source 
of  the  main  geomagnetic  field  is  beneath  the  earth's  crust. 
Approximately  907o  of  the  main  field  exhibits  characteristics 
roughly  equivalent  to  the  field  of  a  terrestrially  centered 
short  bar  magnet  or  dipole  inclined  at  11.5°  to  the  axis  of 
rotation.   See  Figure  2.1. 


Figure  2.1   Dipole  Appearance  of  Geomagnetic  Field 
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The  field  strength  is  approximately  60  kilogammas  in 
the  polar  regions,  30  kilogammas  in  the  equatorial  region, 
and  50  kilogammas  in  the  neighborhood  of  37  North  latitude 
on  the  west  coast  of  the  United  States.   Figure  2.2  shows 
a  world  chart  of  the  total  intensity  of  the  surface  field. 
The  main  field  has  a  slow  secular  change  of  about  .  17o  per 
year  and  it  is  a  regional  phenomenon  which  may  reach  . 3% 
in  certain  regions. 

B.   UNITS 

The  customary  unit  of  measurement  is  the  gamma,  regarded 
interchangeably  as  a  unit  of  induction  (gauss)  or  of  magnetic 
intensity  (oersted) . 


-S        •   -9 
1  gamma  is  10   gauss  (10  v   tesla) 


or 


1  gamma  is  10   oersted  (10/4  it  ampere  turn/meter) 

C .   ELEMENTS 

The  various  elements  of  the  magnetic  field  are  shown 
in  Figure  2.3.   They  are  defined  as: 
X  North  South  Component 


Air  Force  Cambridge  Research  Laboratories,  Handbook  of 
Geophysics  and  Space  Environments,  p,  11-11,  1965. 
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Y  East  West  Component 

Z  Vertical  Intensity 

D  Declination 

H  Intensity  of  the  Horizontal  Component 

F  Total  Intensity 

I  Inclination  or  Dip 


EAST  (-*) 


^  NORTH  (+) 


Figure  2.3   Elements  of  the  Field 


In  terms  of  the  more  familiar  bar  magnet:   A  magnet 
perfectly  free  to  turn  in  any  direction  would  show  the 
vector  direction  of  the  field.   It  is  almost  impossible 
to  suspend  a  magnet  like  that  so  two  are  used.   One  with 


17 


a  vertical  axis  (compass)  and  one  with  a  horizontal  axis 
(dip  needle) .   The  compass  shows  the  direction  of  the 
horizontal  component  H,  and  when  properly  oriented  the  dip 
needle  shows  the  inclination  or  dip. 

D.   TIME  VARIATIONS 

The  time  variations  of  the  geomagnetic  field  can  be 
broadly  categorized  into  quiet  variation  fields,  disturbed 
variation  fields,  and  micropulsations . 

1.   Quiet  Variation  Fields 

Geomagnetic  records  occasionally  show  a  very  smooth 
trace  that  clearly  indicates  a  pattern  of  daily  variation 
with  respect  to  solar  local  time.   This  pattern  is  named 
"solar  quiet  day  variation"  and  denoted  Sq.   There  is  also 
a  lunar  daily  variation,  denoted  L,  but  it  is  much  weaker 
and  more  variable.   Both  the  Sq  and  L  variations  result 
from  the  dynamo  action  of  the  atmosphere,  a  complex  phenom- 
enon of  convective  motions  caused  by  the  solar  heating  of 
the  upper  atmosphere  under  the  influence  of  the  earth's 
magnetic  field.   These  convective  motions  induce  ionospheric 
currents  and  electric  currents  in  the  earth  itself,  both 
of  which  cause  magnetic  variations.   The  effects  of  electric 
current  induction  within  the  earth  are  sometimes  referred 
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to  as  geological  disturbances  because  the  amount  of  induc- 
tion is  dependent  on  the  composition  of  the  earth's  crust. 
For  example,  the  intensity  of  induced  currents  in  the  ocean 
is  higher  than  those  in  the  continental  crustal  layers. 

Figure  2.4  shows  the  Sq  variation  recorded  in 
Monterey  during  a  quiet  period. 

2.   Disturbed  Variation  Fields 

Any  fluctuation  other  than  the  quiet  day  solar  and 
lunar  variations  can  be  called  a  magnetic  disturbance. 
Other  than  disturbances  directly  attributable  to  a  solar 
flare  or  fluctuations  produced  by  upper  atmospheric  irreg- 
ular motions,  all  are  caused  originally  by  disturbed  solar 
plasma  and  may  therefore  be  generally  classified  as  one 
family.   The  term  magnetic  storm  is  reserved  to  apply  to 
a  relatively  severe,  long  lasting  disturbance  with  recog- 
nizable features,  but  all  other  irregular  variations 
deviating  from  the  Sq  and  L  variations  are  termed  simply 
magnetic  disturbances.   They  have  durations  lasting  from 
100  hours  to  several  minutes.   Figure  2.5  shows  a  magnetic 
disturbance  recorded  in  Monterey  that  immediately  followed 
the  Sq  variations  shown  in  Figure  2.4. 
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3.   Geomagnetic  Micropulsations 

Geomagnetic  micropulsations  are  field  variations 
with  periods  ranging  from  .2  seconds  to  approximately  10 
minutes.   A  terminology  recommended  by  the  International 
Association  of  Geomagnetism  and  Aeronomy  is  pc  for  pulsa- 
tion continuous  and  pi  for  pulsation  irregular.   The  period 
range  for  these  are: 

pc  1   .2-5  sec 

pc  2   5-10  sec 

pc  3   10-45  sec 

pc  4   45-150  sec 

pc  5   150-600  sec 

pi  1   1-40  sec 

pi  2   40-150  sec 
A  general  picture  of  the  geomagnetic  field  in  the 
lower  frequencies  is  shown  in  Figure  2.6. 

Pc  1  pulsations  have  very  small  amplitude,  between 
.01  and  .1  gamma  and  have  a  sinusoidal  and  beating  form. 
The  pulsations  appear  in  bursts  and  usually  last  approxi- 
mately 30  minutes.   They  have  been  identified  by  their 
appearance  and  because  there  is  a  direct  relationship 
between  their  reoccurrence  period  and  mid-frequency.   No 
seasonal  or  latitude  dependence  has  been  identified  but 
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Figure  2.6   Geomagnetic  Micropulsation  Spectrum 

there  is  a  diurnal  variation  in  the  mid-frequency  events, 
high  frequency  events  occur  in  post  midnight  hours,  and 
low  frequency  events  in  the  early  afternoon.   Figure  2.7 
is  a  record  of  a  pc  1  micropulsation  observed  in  Monterey. 

Figure  2.8  shows  the  typical  form  of  the  pc  2,3 
pulsations.   They  also  exhibit  a  sinusoidal  and  beating 
character  similar  to  the  pc  1  type  but  their  amplitude  is 
slightly  greater,  usually  .5  gamma.   The  duration  of  the 
signals  seems  to  be  related  only  to  the  sensitivity  of  the 
detector,  which  has  lead  to  the  hypothesis  that  when  no 
signal  is  observed  the  pulsation  is  not  completely  absent 
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Figure  2.7   Pel  Micropulsation; 
1.8  gammas  full  scale;  45  second  record 


Figure  2.8   Pc  2,3  Micropulsation; 
1.8  gammas  full  scale;  45  second  record 
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but  rather  too  weak  to  be  detected.   Pc  2,3  events  occur 
mainly  during  daylight  hours  and  have  a  maximum  near  noon. 
The  average  period  undergoes  a  diurnal  variation,  being 
a  maximum  at  midday.   They  also  have  a  27  day  (solar  rota- 
tion) periodicity  and  a  slight  seasonal  variation. 

As  seen  in  the  typical  pc  4  record  shown  in  Figure 
2.9,  these  micropulsations  slowly  rise  to  a  maximum  of 
between  5  and  20  gammas  and  then  symmetrically  decay.   The 
average  signal  lasts  about  one  hour.   In  mid-latitudes  a 
daytime  occurrence  pattern  can  be  observed  but  the  period 
is  shorter  in  the  daytime  than  at  night,  the  opposite  of 
the  pc  2,3. 

The  longest  period  variations,  pc  5,  display  ampli- 
tudes from  10  to  hundreds  of  gammas .   Their  form  is  again 
sinusoidal  and  they  may  last  from  several  minutes  to  several 
hours.   The  largest  amplitudes  are  observed  in  the  auroral 
zone  and  the  occurrences  have  a  direct  relationship  with 
auroral  zone  phenomena. 

The  class  of  micropulsations  in  the  2  to  .005  Hz 
range  with  irregular  form  and  a  close  relationship  to  upper 
atmospheric  and  auroral  zone  phenomena  have  been  designated 
pi  1,2.   Their  amplitudes  vary  greatly,  the  larger  sometimes 
being  called  micropulsation  storms.   Figure  2.10  shows  the 
irregular  appearance. 
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Figure  2.9   Pc4  Micropulsation; 
7  gammas  full  scale;  6  minute  record 
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Figure  2.10   Pi2  Micropulsation; 
7  gammas  full  scale;  6  minute  record 
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III.   MEASURING  INSTRUMENTS 

Today  a  wide  range  of  devices  are  available  for  measur- 
ing the  magnetic  field,  ranging  from  the  extremely  old 
magnetized  needle  to  the  newer  superconducting  quantum 
types.   The  major  considerations  in  the  selection  of  a 
device  in  a  non-observatory  application  are  sensitivity, 
accuracy,  frequency  response,  mobility,  and  reliability. 
Three  types  of  magnetometers  that  have  developed  with  these 
considerations  as  guidelines  are  the  fluxgate,  optically 
pumped  resonance,  and  superconducting  quantum.   They  will 
be  described  and  a  comparison  of  the  different  varieties 
and  other  types  is  given  in  Table  I. 

A.   FLUXGATE 

A  fluxgate  magnetometer,  sometimes  referred  to  as  a 
saturable  core  or  saturable  inductor,  measures  magnetic 
fields  by  utilizing  the  variation  of  permeability  with 
intensity  of  magnetization  in  a  ferromagnetic  core.   An 
easily  saturable  ferromagnetic  core  of  high  permeability 
is  driven  cyclically  into  saturation  by  an  applied  magnetic 
field.   In  the  absence  of  an  external  field,  which  is 
usually  dc  or  very  low  ac  and  small  with  respect  to  the 
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applied  field,  the  induced  voltage  in  the  sense  windings 
will  be  symmetrical,  i.e.,  contain  only  odd  harmonics  of 
the  fundamental  driving  current  that  produced  the  applied 
field.   In  the  presence  of  an  external  field,  the  induced 
voltage  becomes  asymmetrical,  the  asymmetry  and  the  corre- 
sponding even  harmonics  being  proportional  to  the  amplitude 
of  the  component  of  the  external  field  parallel  to  the  axis 
of  the  core. 

The  fluxgate  magnetometer  has  been  used  widely  in 
various  applications  where  sensitivity,  stability,  and 
absolute  accuracy  has  been  sacrificed  for  the  features  of 
simplicity,  low  cost,  low  power,  and  ruggedness.   Since  it 
is  a  directional  device,  measuring  only  the  component 
parallel  to  the  axis  of  the  core,  a  triaxial  arrangement 
of  three  cores  is  necessary  for  total  field  measurements. 

B.   OPTICALLY  PUMPED  RESONANCE 

The  optically  pumped  resonance  magnetometers  (chiefly 
Rubidium,  Cesium,  and  Helium)  measure  the  total  intensity 
of  the  magnetic  field  by  an  extremely  accurate  detection  of 
the  Zeeman  effect.   Due  to  the  Zeeman  effect  the  energy 
levels  associated  with  a  valence  electron  become  split  into 
sublevels  whose  separations  are  dependent  on  the  total 
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intensity  of  the  magnetic  field.   The  detection  of  this 
splitting  is  based  on  the  principles  of  optical  pumping 
and  monitoring.   Optical  pumping  is  described  in  a  tutorial 
manner  in  Reference  4  and  for  the  Cesium  atom  in  Appendix  A. 
Figure  3.1  shows  a  simplified  block  diagram  of  a  typical 
pumping  apparatus.  A  light  source  producing  an  intense  col- 
limated  beam  of  light  is  first  filtered  and  circularly 
polarized  and  then  passed  through  a  resonance  cell  containing 
the  vapor  to  be  pumped.   As  the  light  is  absorbed  and  emitted 
by  the  sample  atoms,  optical  pumping  proceeds  and  results 
in  a  steady  state  non-equilibrium  population  distribution 
over  the  magnetic  states  of  one  of  the  ground  levels.   In 
this  pumped  condition  the  sample  is  prepared  for  resonance 
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Figure  3.1   Typical  Pumping  Apparatus 
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and  is  characterized  by  a  net  magnetization  M  =  Mz  along 
the  ambient  magnetic  field  direction.   Resonance  transitions 
between  the  Zeeman  states  of  the  pumped  level  are  induced 
by  subjecting  the  sample  to  an  RF  field  of  frequency  f-i 
sufficiently  close  to  the  resonant  frequency  £q   and  with 
component  H-j_  perpendicular  to  the  ambient  field  Hq  =  Hz. 
The  RF  drive  field  introduces  Zeeman  coherence,  and  accord- 
ing to  Bloch's  equations  the  steady  state  behavior  of  the 
magnetization  M  changes  and  is  no  longer  parallel  to  Hq  but 
now  has  a  transverse  component  M.  precessing  around  Hq  at 
the  frequency  of  H,  .   Consequently  the  M2  component  along 
Hq  decreases,  and  the  light  transmission  also.   This  varia- 
tion of  light  intensity  with  the  applied  RF,  being  a  maximum 
at  resonance,  can  be  used  as  a  detection  method,  or  the 
photocell  can  be  made  to  observe  the  signal  corresponding 
to  the  M  component,  which  is  precessing  or  rotating  at  the 
same  RF  applied  which  induced  the  resonance.   The  photocell 
output  is  then  at  the  resonant  frequency.   Appendix  B 
details  the  two  detection  systems  commonly  used.   The 
resonant  frequency  is  proportional  through  a  constant  to 
the  magnetic  field  (4.66  Hz/gamma  for  Rb ;  3.499  Hz/gamma 
for  Cs;  28  Hz/gamma  for  He)  and  thus  the  field  measurement. 
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Figures  2.4  through  2.10  were  all  recorded  using  an 
optically  pumped  magnetometer  and  the  appearance  of  the 
micropulsations  differ  somewhat  from  those  one  would  find 
in  a  reference  handbook.   This  illustrates  an  undesirable 
feature  of  this  magnetometer  when  used  in  micropulsation 
work.   It  is  equally  sensitive  to  all  frequencies  and  the 
larger  amplitude,  lower  frequency  changes  such  as  the  di- 
urnal variation  continually  drive  the  highly  sensitive 
instrument  off scale.   The  result  is  a  sloping  and  discon- 
tinuous output  that  is  difficult  to  interpret. 

C.   CRYOGENIC  SQUID 

A  cryogenic  SQUID  (Superconducting  Quantum  Interference 
Device)  magnetometer  consists  of  three  basic  parts:   1)  a 
superconducting  ring  embodying  a  Josephson  junction;  2)  an 
externally  generated  HF  field  used  to  provide  a  small  field 
variation  superimposed  on  the  ambient  dc  field  and  therefore 
to  control  the  amount  of  flux  through  the  ring;  3)  an  elec- 
tronic detection  system  to  observe  the  impedance  of  the 
coil  and  make  the  output  response  a  linear  function  of  the 
ambient  field. 

The  operation  utilizes  the  principle  that  if  the  flux 
through  the  ring  is  varied  on  either  side  of  a  minimum  or 
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maximum  energy  value,  associated  with  the  Josephson  coupling 
energy,  the  SQUID  will  generate  a  field  which  will  respec- 
tively oppose  or  add  to  the  variation.   The  energy  contribu- 
tion of  the  SQUID  aiding  or  opposing  the  applied  field 
variation  according  to  whether  the  ambient  dc  field  amounts 
to  a  minimum  or  maximum  energy  value  constitutes  a  para- 
metric control  of  the  inductance  measured  at  the  coil 
terminals.   In  other  words,  the  SQUID  alternately  behaves 
diamagnetically  and  paramagnetically  and  will  decrease  and 
increase  the  inductive  reactance. 

Through  the  use  of  a  detector  to  monitor  the  coil  imped- 
ance one  can  observe  the  variation  of  impedance  with  ambient 
magnetic  field.   Because  the  impedance  variation  is  periodic 
a  phase  locked  feedback  loop  is  normally  included  to  obtain 
a  response  linearly  related  to  the  field. 

Because  SQUID  compares  the  energy  values  of  the  ambient 
and  weak  superimposed  fields,  it  measures  only  a  change  in 
field,  not  the  total  intensity. 
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D.   GRADIOMETER 

A  gradiometer,  in  theory,  would  measure  the  space 
gradients  of  the  magnetic  field,  that  is: 


<3x  ,  3x_  ,  dx 

8x  dy  ,       dz 

3y  ,  3y  ,    dY_ 

3x  8y  8z 

dz_  ,  8z_  ,  9z 

8x  8y  3z 


Only  five  of  these  are  independent  because  X,  Y,  and  Z 
are  the  derivatives  of  a  potential  and  because  of  Laplace's 
equation.   In  practice  gradiometers  are  constructed  by 
taking  the  finite  difference  of  two  magnetometers,  there- 
fore cancelling  the  ambient  field  that  is  approximately 
equal  at  both  sensors.   A  local  change  or  anomaly  that  is 
closer  to  one  of  the  two  sensors  is  not  cancelled  and 
therefore  recorded. 

The  disadvantage  of  the  gradiometer  is  that  with  respect 
to  a  single  sensor  magnetometer  most  magnetic  anomalies  fall 
off  as  the  inverse  cube,  1/r  where  r  is  the  distance  from 
the  anomaly  to  the  sensor,  but  the  gradiometer  anomaly  falls 
off  as  the  inverse  fourth  power,  1/r  where  r  is  the 
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distance  from  the  anomaly  to  the  center  of  the  sensor 
configuration . 

In  Reference  10,  Langan,  Lynch,  and  Slack  reported  a 
gradiometer  used  in  geophysical  prospecting  that  measured 
the  vertical  gradient  A  H/az  where  H  was  the  magnitude  of 
the  anomaly  and  z  the  distance  to  the  anomaly.   It  was 
constructed  of  two  optically  pumped  sensors  vertically 
separated  and  towed  by  a  helicopter. 

Both  the  optically  pumped  and  SQUID  magnetometers  are 
easily  converted  to  gradiometers  because  of  the  form  of 
their  output,  but  they  will  measure  different  gradients. 
The  optically  pumped  gradiometer,  constructed  of  total 
field  magnetometer  sensors  will  measure  the  total  field 
gradient   AF/Az.   A  SQUID  gradiometer,  constructed  of 
devices  sensitive  only  to  a  field  change  will  measure  a 
field  component  gradient  AZ/  Az,  for  example. 
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IV.   SYSTEM  DESCRIPTION 

The  overall  system  can  be  broken  into  three  subsystems: 
1)  the  data  collection  subsystem;  2)  the  periodmeter  and 
interface  subsystem;  3)  the  processing  and  display  subsystem, 


A.   DATA  COLLECTION  SUBSYSTEM 

The  magnetometer  used  was  the  Varian  V-4938  Cesium 
system.   It  used  an  optically  pumped  cesium  sensor  of  the 
auto-oscillating  type.   A  block  diagram  of  the  Varian  sys- 
tem is  shown  in  Figure  4.1. 
SENSOR 


200  FT 


LARMOR 

FREQUENCY 
-176  kHz 


READOUT  UNIT 


dc 
VOLT 


PEN 
RE- 
CORDER 


Figure  4.1   Varian  Magnetometer 
The  sensor  head  produces  the  RF  signal  (Larmor  frequency) 
directly  proportional  to  the  total  magnetic  field  intensity. 
This  signal  together  with  the  input  power  for  the  sensor  is 
conducted  along  a  200  foot  coaxial  cable  to  the  main  readout 
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unit.   The  readout  unit  converts  the  Larmor  frequency  to  a 
highly  stable  dc  voltage  which  drives  the  pen  recorder. 
The  system  was  located  in  a  wooded  area  approximately  one 
mile  from  the  Naval  Postgraduate  School,  which  provided  a 
near  noise  free  environment. 

A  more  complete  block  diagram  of  the  V-4938  system  is 
shown  in  Figure  4.2.   The  RF  Larmor  frequency,  which  was 
approximately  176  kHz,  from  the  sensor  is  mixed  in  the  read- 
out unit  to  the  100  Hz  to  3  kHz  range  and  then  applied  to 
a  highly  stable  frequency  to  voltage  converter.   This  dc 
voltage,  which  now  represents  the  field  value,  is  passed 
through  biasing  circuitry  to  interface  it  to  the  pen  recorder. 

The  system  configuration  used  for  data  collection  differed 
from  the  Varian  recommended  system  of  Figure  4.2.   This  was 
necessary  because  of  the  requirement  to  record  the  output 
on  magnetic  tape.   The  actual  system  used  is  shown  in  Figure 
4.3.   This  configuration  was  used  because  the  mixer  output 
of  the  readout  unit  was  a  square  wave  and  its  frequency 
content  exceeded  the  5  KHz  range  of  the  tape  recorder.   The 
coupler  unit  shown  connects  between  the  readout  unit  and 
the  sensor  to  provide  dc  power  to  the  sensor  and  it  also 
contains  a  mixer,  but  its  output  is  a  sinusoid  rather  than 
a  square  wave  and  recorded  nicely.   The  coupler  units  primary 
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Figure  4.2   Magnetometer  Signal  Flow 
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Figure  4.3   Magnetometer  Data  Collection  Configuration 
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purpose  is  to  convert  the  system  to  a  gradiometer  which  is 
shown  in  Figure  4.4. 


BIAS 


SENSOR  #1 


COUPLER 

#1 


MIXER  OUTPUT 
GRADIOMETER 


READOUT 


FREQ. 

TO 
VOLT 


B 


SENSOR  #2 


COUPLER  #2 


1  MIXER  OUTPUT 
MAGNETOMETER 


PEN 
RE- 
CORDER 
(GRADIO 
METER) 


Figure  4.4   Gradiometer  Data  Collection  Configuration 

Again  this  is  not  the  configuration  recommended  but 
this  provided  both  a  magnetometer  and  gradiometer  output 
suitable  for  recording  and  still  allowed  the  pen  recorder 
to  monitor  the  gradiometer  signal.   The  bias  applied  to 
the  one  sensor  was  required  to  generate  a  difference 
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frequency  between  the  sensors  sufficient  to  drive  the  mixer 
circuitry.   The  sensors  were  located  approximately  180  ft 
apart  resulting  in  a  Larmor  frequency  difference  of  only 
approximately  30  Hz,  and  the  mixers  required  a  difference 
of  at  least  100  Hz  and  optimally  between  1  and  3  kHz. 
Figure  4.5  shows  the  readout  units  mixer  output  with  an 
insufficient  difference  frequency  and  Figure  4.6  shows  the 
coupler  units  sinusoidal  mixer  output. 

One  other  system  characteristic  imposed  a  limitation, 
that  was  the  orientation  sensitivity  of  the  single  cell 
sensor.   For  optimum  performance  the  sensor  head  needed  to 
be  oriented  both  toward  magnetic  North  and  45°  to  the  total 
field  vector.   Figure  4.7  shows  the  non-metallic  housing 
and  Figure  4.8  the  geometry. 

Figure  4.9  shows  the  system  components  and  equipment 
used  to  monitor  and  record  the  data. 
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Figure  4.5   Distorted  Readout  Unit  Mixer  Output 


Figure  4.6   Coupler  Unit  Mixer  Output 
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Figure  4.7 


Non-Metallic  Housing  Used  to  Fulfill 
Orientation  Requirements 


42 


FIELD 
VECTOR 


NORTH 


„-.  3PTICAL 
AXIS 


Figure  4.8   Sensor  Housing  Geometry 
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B.   PERIODMETER  AND  INTERFACE  SUBSYSTEM 

The  frequency  response  of  the  magnetometer  when  used 
with  the  readout  unit  is  limited  to  20  Hz  by  the  discrimin- 
ator and  even  more  by  the  pen  recorder  response.   The 
theoretical  frequency  response  of  the  sensor  is  limited 
only  by  the  bandwidths  of  the  amplifier  and  feedback  loop, 
which  are  extremely  wide.   To  take  advantage  of  this  the 
mixer  output  was  recorded  on  an  HP-3960  FM  tape  recorder. 
This  signal  still  contained  the  high  frequency  components 
present  in  the  Larmor  signal.   A  periodmeter  and  interface 
unit  were  constructed  to  measure  the  period  of  successive 
cycles  and  interface  this  digital  count  into  the  XDS-9300 
computer,  for  further  display  on  the  AGT-10  graphics  console. 

Figure  4.10  is  a  block  diagram  of  the  periodmeter.   The 
data  signal  is  passed  through  a  simple  RC  high  pass  filter 
to  reject  any  60  Hz  modulation  that  might  be  present  and 
then  input  to  a  zero  crossing  detector.   The  falling  edge 
of  this  is  used  to  trigger  a  one-shot  that  inhibits  the 
otherwise  free  running  4  MHz  count.   The  one-shot  pulse 
is  inverted  and  delayed  and  the  rising  edge  used  to  trigger 
flip  flop  type  memories.   Once  the  count  is  latched  another 
one  shot  is  used  to  clear  the  count  and  provide  a  signal 
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indicating  the  cycles  completion.   The  cycle  then  repeats. 
The  actual  circuits  and  timing  diagrams  are  presented  in 
Appendix  C. 

A  direct  access  channel  into  the  XDS-9300  was  used  to 
input  the  binary  count  directly  into  sequential  memory  loca- 
tions.  Figure  4.11  is  a  block  diagram  of  the  periodmeter 
and  interface.   To  implement  the  input  the  following  sequence 
of  instructions  were  used: 

EOM  032400 

SKS   032400 

BRU  $-1 

PIN  ITEM 
With  the  execution  of  the  EOM  instruction  the  contents  of 
the  control  register  are  available  at  the  output  connector. 
This  is  decoded  and  toggles  a  flip  flop  to  alert  the  period- 
meter  that  the  XDS  is  ready  to  accept  an  input.   This  EOM 
alert  and  the  periodmeter ' s  completion  signal  are  logically 
combined  and  send  a  ready  signal  to  the  XDS-9300.   Upon 
receipt  of  this  ready  signal  and  the  execution  of  the  PIN 
the  data  is  transferred  and  stored  in  the  location  specified 
in  the  PIN  instruction.   The  SKS  instruction  is  used  to  sense 
if  both  the  computer  and  periodmeter  are  ready.   If  both  the 
EOM  alert  and  completion  signal  are  not  present,  the  next 
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instruction,  branch  unconditionally  back  to  the  SKS ,  is 
executed.   This  looping  continues  until  both  systems  are 
ready.   If  both  are  ready  the  branch  instruction  is  skipped, 
The  ready  signal  is  also  used  to  gate  the  data  through  in- 
verters because  it  is  inverted  prior  to  being  transferred 
into  the  XDS-9300.   The  actual  circuit  diagrams  are  pre- 
sented in  Appendix  C. 

C.   DISPLAY  SUBSYSTEM 

A  block  diagram  of  the  data  flow  in  the  interface  and 
display  subsystems  is  shown  in  Figure  4.12. 
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Figure  4.12   Data  Flow  From  Periodmeter  Through  Display 

The  physical  block  length  of  data  input  to  the  XDS-9300 
was  12,288  points.   A  mean  and  standard  deviation  of  these 
points  was  found  and  values  not  within  five  standard 
deviations  were  discarded.   This  eliminated  erroneous  data 
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points  caused  by  noise  spikes  either  in  the  computer, 
periodmeter,  or  on  the  recording.   The  block  was  then 
buffered  out  in  groups  of  4096  points  and  written  on 
magnetic  tape.   To  be  displayed  the  data  was  buffered 
back  in,  again  in  groups  of  4096  points,  normalized,  and 
displayed  on  the  AGT-10  graphics  console.   Both  programs 
are  self-explanatory  and  are  presented  in  Appendix  D. 
Signal  flow  diagrams  for  the  two  are  shown  in  Figure  4.13 
and  4 . 14 . 
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Figure  4.13   Signal  Flow  Diagram  of  Input  Program 
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Figure  4.14   Signal  Flow  Diagram  of  Display  Program 
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V.   RESULTS  AND  CONCLUSIONS 

The  system  was  tested  in  three  phases:   1)  laboratory 
signals;  2)  tape  recorded  laboratory  signals;  3)  magnetom- 
eter signals. 

First,  constant  frequency  signals  from  a  laboratory 
generator  were  applied.   The  frequencies  sampled  and  the 
resulting  period  measurement  are  listed  in  Table  2. 

FREQUENCY  Hz    MEAN  PERIOD ^sec   THEORETICAL  PERIOD  ^sec 

1000.0  999.25  1000.00 

1000.1  999.00  999.00 
1000.5  998.75  999.50 
1001.0  998.25  999.00 
1003.0  996.25  997.01 
1005.0  994.25  995.02 

Table  2 


The  .75  usee  difference  between  theoretical  and  actual  is 
due  to  the  clear  and  latch  time  of  the  periodmeter.   These 
measurements  were  within  +1  count  or  +25  usee  with  a  4MHz 
time  base.   This  established  the  resolution  of  the  device 
as  .1  Hz  or  better.   The  outputs  displayed  for  the  sampled 
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frequencies  of  1000.0  Hz,  1000.1  Hz,  and  1001.0  Hz  are 
shown  in  Figures  5.1  -  5.3  and  a  distinguishable  difference 
is  noticeable.   To  insure  the  system  would  follow  a  rapidly 
changing  frequency  a  swept  frequency  signal  was  applied  and 
the  display  shown  in  Figure  5.4. 

The  second  phase  of  testing  was  designed  to  determine 
the  magnitude  of  the  tape  recorder  flutter  and  noise.   It 
involved  recording  the  same  frequencies  used  in  phase  1  and 
comparing  the  results.   The  frequencies,  mean  period  measure- 
ment, and  standard  deviation  for  each  are  listed  in  Table  3. 

FREQUENCY  Hz    MEAN  PERIOD  usee    STANDARD  DEVIATION  usee 

1000.0  999.25  .575 

1000.1  999.00  .700 
1000.5  998.75  .850 
1001.0  998.25  .550 
1003.0  996.25  .650 
1005.0  994.25  .675 

The  mean  period,  calculated  for  12,288  periods,  remained 
the  same.   The  average  standard  deviation  was  .667  usee. 
This  established  the  RMS  noise  level  of  the  system  when 
operated  in  the  1200  Hz  range  to  be  approximately  1  Hz  or 


54 


1000. 


PERIOD 
fJlSEC 


990. 


0 


1.2 


TIME,  SEC 


Figure  5.1  Output  Display  of  Direct  Input 
1000.0  Hz  Laboratory  Signal 


1000. 


PERIOD 
pSEC 


990. 


TIME,  SEC 


Figure  5.2  Output  Display  of  Direct  Input 
1000.1  Hz  Laboratory  Signal 
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Figure  5.3  Output  Display  of  Direct  Input 
1001.0  Hz  Laboratory  Signal 
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Figure  5.4   Output  Display  of  Direct  Swept  Frequency 
Laboratory  Signal 
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.27  gamma.   The  AGT-10  displays  for  the  sampled  frequencies 
of  1000.0  Hz  and  1000.1  Hz  are  shown  in  Figures  5.5  and  5.6. 
The  previously  discernible  difference  is  now  buried  in  tape 
recorder  noise.   Figure  5.7  shows  the  frequency  transition 
from  1001.0  to  1003.0  Hz  and  this  2  Hz  difference  is  easily 
dis  tinguishable . 

The  final  phase  of  testing  involved  introducing  a  known 
magnetic  disturbance  and  analyzing  the  magnetometer  output. 
A  72  turn  Helmholtz  coil,  in  the  configuration  of  Figure  5.8, 
was  used  to  introduce  the  disturbance.   The  sensor  was 
placed  on  the  axis  of  the  coil  and  the  coil  aligned  with 
the  total  field  vector.   This  insured  the  induced  field  of 
the  coil  would  directly  aid  or  oppose  the  ambient  field 
and  allowed  changes  of  predetermined  amounts .   With  the 
simplified  geometry 

H  =  N  I  R2    4tt 
2  r3     10_z 

where  N  =  number  of  turns  in  the  coil 

I  =  current  in  the  coil 

R  =  radius  of  the  coil 

r  =  distance  from  coil  to  sensor 

<T     was  required  for  conversion  from 
10     ampere-turns  per  meter  to  gammas. 
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Output  Display  of  Recorded 
1000.0  Hz  Laboratory  Signal 
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Figure  5.6  Output  Display  of  Recorded 
1000.1  Hz  Laboratory  Signal 
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Figure  5.7  Output  Display  of  RecordedlOOl.O 
and  1003.0  Hz  Laboratory  Signal 
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Coil  and  Sensor  Geometry  Used 
to  Disturb  the  Field 


60 


The  coil  utilized  had  a  mean  radius  of  33  cm  and  was  placed 
3.5  m  from  the  sensor.   Table  4  tabulates  the  various 
currents  and  observed  versus  predicted  changes.   Figures 
5.9  through  5.30  show  the  magnetometer  pen  recorder  output 
and  a  selected  time  frame  of  the  corresponding  AGT-10  display, 
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Figure  5.9  Pen  Recorder  Output  for  AH=  45  gamma 
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Figure  5.10  AGT-10  Display  for  AH=  45  gamma 
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Figure  5.11  Pen  Recorder  Output  for  AH=  35  gamma 
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Figure  5.12  AGT-10  Display  for  A  H=35  gamma 
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Figure  5.13   Pen  Recorder  Output  for  A  H=  28  gamma 
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Figure  5.14  AGT-10  Display  for  A  H=  28  gamma 
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Figure  5.15   Pen  Recorder  Output  for  AH=  18  gamma 


1000. 


o 
w 
w 


Q 
O 
M 

W 
Pk 


820. 


TIME,  SEC 
Figure  5.16  AGT-10  Display  for  AH=  18  gamma 
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Figure  5.17  Pen  Recorder  Output  for  AH=  11  gamma 
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Figure  5.18  AGT-10  Display  for  A  H=  11  gamma 
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Figure  5.19   Pen  Recorder  Output  for  A  H=  8  gamma 
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Figure  5.20  AGT-10  Display  for  A  H=  8  gamma 
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Figure  5.21   Pen  Recorder  Output  for  A  H=  4  gamma 
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Figure  5.22  AGT-10  Display  for  AH=  4  gamma 
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Figure  5.23  Pen  Recorder  Output  for  AH=  2  gamma 
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Figure  5.24  AGT-10  Display  for  A H=  2  gamma 
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Figure  5.25   Pen  Recorder  Output  for  AH=  1  gamma 
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Figure  5.26  AGT  Display  for  AH=  1  gamma 
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Figure  5.27  Pen  Recorder  Output  for  AH=  .6  gamma 
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Figure  5.28  AGT-10  Display  for  AH=  .6  gamma 
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To  obtain  a  precise  measurement  of  the  magnetometer's 
response  time  the  AGT-10  display  was  expanded.   Figure  5.29 
shows  the  rising  edge  and  Figure  5.30  the  falling  edge.   As 
seen  the  rising  edge  required  approximately  three  periods 
or  2.5  msec  to  respond  and  the  falling  edge  only  one  period 
or  .8  msec. 
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Figure  5.29   Expanded  View  of  Rising  Edge 
of  Magnetometer  Response 


72 


1000. 


w 


Q 
O 

M 

w 

p^ 


820. 


0. 


TIME,  SEC 


1.2 


Figure  5.30   Expanded  View  of  Falling  Edge 
of  Magnetometer  Response 
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The  periodmeter  and  interface  unit  measured  the  period 
of  the  signal  with  an  absolute  accuracy  of  .75  usee  and 
properly  input  the  value  to  the  computer.   The  resolution 
of  the  device  was  .25  usee  using  a  4  MHz  time  base  and 
this  could  be  increased  by  using  a  faster  time  base.   The 
system  was  successfully  operated  with  a  10  MHz  clock  fre- 
quency and  resolution  of  .1  usee.   The  use  of  an  increased 
clock  frequency  or  averaging  of  periods  to  increase  the 
periodmeters  resolution  would  not  improve  the  system  sensi- 
tivity.  The  overall  system  sensitivity  is  limited  by  the 
flutter  and  noise  of  the  tape  recorder.   Two  methods  were 
attempted  to  overcome  this:   1)  the  use  of  a  synchronized 
oscillator  constructed  of  two  multivibrators  synchronized 
with  the  zero  crossings  of  the  input  signal;  2)  recording 
a  1  kHz  reference  signal  on  a  second  channel  of  the 
recorder  and  multiplying  it  up  for  use  as  a  time  base. 
The  pulse  width  instability  of  the  multivibrators  was 
greater  than  the  flutter  noise  so  the  synchronized  oscillator 
was  unsatisfactory.   A  phase  lock  loop  type  frequency  multi- 
plier was  constructed  to  multiply  the  1  kHz  to  1  MHz  but 
the  delay  time  through  the  loop  was  too  great  to  cancel  the 
flutter  for  a  period  to  period  measurement. 

To  achieve  greater  system  sensitivity  and  resolution  some 
means  of  flutter  cancellation  is  necessary. 
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APPENDIX  A 
OPTICAL  PUMPING  OF  THE  CESIUM  ATOM 

Figure  A-l  shows  the  energy  levels  of  the  cesium  atom. 

Absorption  of  photons  of  light  by  the  cesium  valence  elec- 

o 
trons  cause  their  energy  to  be  increased  from  the  Sl  ground 

o 
state  to  the  Pi  state.   Due  to  the  light's  polarization 

the  magnitude  of  the  quantum  number  m  is  also  increased  by 

one.   An  excited  atom  quickly  emits  its  energy  and  returns 

to  any  of  the  ground  state  sublevels.   The  absorption  and 

re-emission  is  repeated  resulting  in  a  shift  to  the  m=+4 

2 
sublevel.   Since  the  Pi  state  contains  no  m  level  above 

% 

2 
+4,  transitions  out  of  the  S^,  m=+4  level  are  not  possible. 

This  is  precisely  optical  pumping  and  under  these  conditions 

there  can  be  no  further  absorption  of  light. 
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Figure  A-l   Energy  Levels  of  Cesium  Atom 
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APPENDIX  B 
OPTICAL  DETECTION  SYSTEMS 

The  first  method  used  to  obtain  a  continuous  record  of 
the  total  field  intensity  used  a  control  loop  arrangement 
and  was  therefore  called  a  Controlled  Magnetometer.   A 
block  diagram  of  it  is  shown  in  Figure  B-l. 
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Figure  B-l  Controlled  Magnetometer 

A  small  frequency  modulation  is  applied  to  the  local  RF 
oscillator  to  sweep  it  back  and  forth  through  resonance. 
Owing  to  this  sweeping  the  signal  collected  at  the  photo- 
cell is  also  modulated  at  the  same  frequency,  but  its  phase 
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with  respect  to  the  initial  modulating  signal  depends  on 
the  position  of  the  local  RF  with  respect  to  the  resonant 
frequency.   The  photocell  output  and  modulating  signal  are 
phase  detected  and  the  output  is  used  as  an  error  signal 
in  a  control  loop  to  maintain  the  local  RF  oscillator  coinci- 
cent  with  the  resonant  frequency.   This  system  does  not  place 
severe  restrictions  on  the  performance  of  the  photocell  or 
its  amplifier,  has  high  accuracy,  good  stability,  and  good 
low  field  capability. 

One  attribute  the  controlled  system  does  not  have  is  a 
good  response  to  transients.   Figure  B-2  shows  the  block 
diagram  of  a  system  that  has  the  excellent  response  to 
transients  necessary  to  detect  and  observe  micropulsations. 
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Figure  B-2   Self-Oscillating  Magnetometer 
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In  theory  a  second  light  beam  is  placed  at  right  angle  with 
respect  to  the  first  and  the  variation  of  the  M  component 
of  the  magnetization  is  observed.   Since  this  transverse 
component  oscillates  at  the  frequency  of  the  signal  injected 
and  the  amplitude  is  a  maximum  at  resonance,  the  photocell 
is  simply  amplified,  phase  shifted,  and  re-injected  to  the 
coils.   The  system  will  then  oscillate  at  the  resonant 
frequency.   In  practice  it  is  possible  to  use  a  single  light 
beam  inclined  at  45°  to  the  direction  of  the  field  to  be 
measured.   This  places  an  orientation  restriction  on  the 
sensor,  but  systems  employing  two  and  three  resonance  cells 
have  been  developed  that  overcome  this.   This  system  places 
severe  restrictions  on  the  amplifier  and  especially  on  the 
accuracy  of  the  phase  shift,  which  is  required  to  maintain 
oscillation.   It  has  a  low  absolute  accuracy  and  poor  long 
term  stability  but  as  previously  mentioned,  excellent 
response  to  transients. 
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APPENDIX  C 
PKRTODMETER  AND  INTERFACE  CIRCUITS 
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Figure  C-2   Input  Circuitry  Timing  Diagram 
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APPENDIX  D   COMPUTER  PROGRAMS 


-JOB 

-FCRT^AN    LStGCS 

CINENSICN    CGUNK13C00)  ,IVAR<6) 

IMEGER    CCLNT 

CJTA    NSAMP, MULT/4096, 3/ 

LELK=NSAMF*MIJLT 
106    FCPMATC  IX, 12(110)  ) 
9^9    CCNTINLE 
C 

C    INITIALIZE  VARIALBLSS 
C 

IM>=1 

IMN=1GC000 

SLMO. 

StMSC=C. 

C=C. 

IEAC=0 

E*C=C. 
C 

C   INFL1  CATA  FROM  PER  ICCNETER/  INTER  F/!CE 
C 

CC  2  IaltLBLK 
5      ECN  C32400 
S      SKS  C3240C 
S      EFL  $-1 
S      FIN  ITEN 

ICCLNT  =0 

ICCLNT=LANC( ITEM, 177B) 

ICCLNT=IICR{IC0UNT,LANC(LRS(ITEN,1) ,37600  6)  ) 

ICCLNT=LICR(ICOUNT,LANC(LRS(ITEV  ,1),140CCCB)) 

CCLNTd  )=ICCUNT 

2  CCNTINLE 

r 

C  FCFf*  *N  INITIAL  AND  THEN  PRECICE  MEAN  £  SCEV 
C 

CC  3  I=1,LELK 

SlN=SLM+CCLNTd) 

SLMSCsSUMSC+COUNT (I )**2 

C=C+1.0 

3  CCNTINLE 
MEAN=SLN/C 

SCEV=SQRTf ABS(C*SL*SC-SLM**2) )/C 
WFI7£(6,1G0)MEAN,SCEV 
ICC  FCFNAT(*0  MEAN  = » , F12. 4/ 
1*  SCEV=',F12.4) 
CC  4  I=1,LBLK 

IF(ABS<CCUNT(I)-MEAN).LE.5.*S0EV)    GC    TO    4 
SIMSUM-CCONT(I) 
SINSG^SL NSC-COUNT (I )**2 
E/C»EAC4l. 
A    CCNTINUE 
C=C-BAC 
MAK*SU*»/C 

SCEV=SQFTIAES(C*SLMSC-SLM**2)  )/C 
fcFITE(6,101)MEAN,SCEV 

101  FCFMATCO    TRUE    MEAN= • , F 12.4/ 
1»    TRUE    SDEV=»,F12.4) 

C  FCP*<    A    NEk    ARRAY    CONSISTING    OF    POINTS    WITHIN     10    SDEV 

C 

I  E AC=0 

N£U=0 

CC    t    I=1,LBLK 

IF (ABS(CCLNT(I ) -MEAN ) .L E. 1Q.*SDE V )NEW=NEkI+ 1 ; GC  TC  7 

IEAC=IEAC*1 

VsFI7E(6,lC2)  IBAD, I, CCLNTd  ) 

102  FCFfATCO  THIS  IS  BAC  PCINT  N0»,I5/ 
1'  ITS  OLD  SUBSCRIPT  WAS1 f  18/ 

2*  ITS  VALUE  WA.S',110) 
GC  TC  6 
7  CCINT(NEWI  )=COUNTU) 
t  CCNT1NUE 
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CC  15  I=NEWI+1,LBLK 
CCLNTCI  )  =  C 
If  CCNTINUE 

hF!TE(6fl03)  NEW! 
103  FCFNAT( «C  NEWI*« ,110) 

C  ZEFC  FILL  AFRAY  FROM  LAST  PCINT  TC  ENC  OF  BLCCK 
C 


c 
c 
c 


c 
c 
c 


CC    6    I=NEWI+1,LBLK 
CCLMC  )=0 
E    CCNTIME 

FINC    A    MINIMIS    AND    MAXIMUN    VALUE 

CC    9    1*1, NEW! 
I*IN=AMINO(CGUNT(I>fIHIN) 
IMX=AMAXQ(COUNTm,IMAX) 
9    CONTINUE 

VvPITE  (6  ,1G4)IMIN,!NAX 
10^    FCFNATC'G    MINIMUM*  •  ,1 10/ 
1»    MXIMLM*1  ,110) 

ELFFEFCUT    TC    TAPE    THE    ARRAY    IVAR    CCNTAINING    VARIABLES 

XWRf  1)*NEV>X 

I\AF(2)=MEAN 

IVAR(3)=SCEV 

IWRU)*IMXN 

IV4P(5)«IPAX 

WFITE(6,10S)(IVAR<I  )  ,1 =1, 5 ) 
1CS    FCRNATt  1X.5I10) 

CLTFLT1102)  •WRITE    CN    TAPE1 

REAC(1C1,601)IANSS 
601    FCFNAKA4) 

IFCIANSS.NE.3HY5S)    GO    TC    60  3 

C4LL    BUFFERCUT(1,1*IVAR,5,IND) 

11  IFUNC.EQ.l)    GO    TC    11 
C 

C   ELFFEFCUT  TC  TAPE  THE  ARRAY  CCLNT  IN  GROUPS  CF  4096 
C 

CC    12    I»1,NULT 

J  =  I-1 

M«J*KSAf<F)  +  l 

CALL    BUFFEROUTd,  1,  COUNT*  K)  ,NSAPF,IND) 

12  1FUNC.EC.1)    GO    TC     12 
12    CONTINUE 

6C3    CCN7INUE 

CLTFUTU02)     'READ    ANOTHER    14336' 
RE*CUG1,600)IANS 
600    FCFNAT(A4) 

IF(IANS.EC.3HYES)    GO    TO    999 
EI^C 
-^LC4C    >F,NAP 
^CAT* 
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PATCH 
>>C*7 
07067 

ENC 
AGT 
JCE 
FCFTP 


CC1C6711 


AN    LS,GC 

C1NENSICN    I  AX  IS ( 10 )» INC  AT  A (4105) ,CATA(41Q5 

ICCIR(5),ITCIR(15),ITEXT(24),T(2) 

INTEGER    RSSP(4) 

CATA    NSAiMP,MULT,NCIS/4G96,3,16/ 

ICEV=1 

LELK=NSAMP*VLL7 

CIV=FLG4T(MCIS) 

NELM1 

XAXIStl  J=IHEAD<0,  10) 

1 1* IS (2)= I  PACK (-.5,. 75,0) 

I/!>IS(2  >=IPACK(-.5,1.,0) 

IA>IS(4)=IPACK<-.5,0.,1) 

I  AX  IS (5 )*IPACK(0.tQ.f 1) 

,-.01,1) 

,0.,GJ 

tC.,1) 

,-.01,1 


)  ,IVAR(6) 


I*XIS(6  )  =  IPACK<0. 
IAXIS(7)=IPACK(0. 
I4XIS(8)=IPACK( .5 
IA>IS<9)=IPACK(.5 
99S  CCNTINLE 
TIfE=C. 
LFEC=NSAMP 
2  C/LL  FNS( ICEV,ISW 
IFUER.NE.G)  CUTPL 
1F(  ISW.NE.Q)  GG  TC 


) 


,IER) 
"(lODIERt 
2 


Fl 


C  BLFFEFIN  ARRAY  IVAR  WITH  VARIABLES 


C/LL    EUFFEPIN(1,1, IV£P 

IF(INC.EC.l)    GO    TC    3 

NEH=1VAR(1) 

N£/N=FLCAT< JVAR(2)  )/4. 

2CEV=FLCAT(XVAR(2 

NIMFL0AT(IVAR(4) 

M /X=FLC/TdVAR(5) 


,5, 


INC) 


)  )/4. 

)/4. 

)/4. 

T  _ 


WFITE(6,8GC)( IVAR (I),  I =1,5) 
800  FCFNATdX,  *NQ.  CF  POINTS* •,  110/ 


1  • 

2  • 
-  i 


N£AN=«,I1G/ 
SCEV=»,I1G/ 
MNIMLM  =  «  ,110/ 


4«  NAXIMLMa'f 1101 
WRITE  MEANfSCEVfMIN»NAX  CN  SCREEN 


2CC 


01 


302 


3C4 


2Qt 


C/ 
IF 
C/ 
IF 

en 

FC 
CA 
IF 

EN 
FC 
CA 
IF 
EN 
FC 
CA 
IF 
EN 
FC 
CA 
IF 
EN 
FC 
14X 


LL  CG 
(IER. 
LL  DT 
(IER. 
CCDE( 
FNAT( 
LL  TE 
(IER. 
CCCE( 
FNAT( 
LL  TE 
(IER. 
CCCE( 
FNAT( 
LL  TE 
(IER. 
CCDEt 
FNAT( 
LL  TE 
(IER. 
CCCE( 
FfAT( 
,'PAX 


INITd 
NE.OJQ 
INITd 
NE.C) 
96,200 
•PE°IO 
XTCdD 
NE.O)0 
96,201 
•T  SEC 
XTCdD 
NE.OC 
96,303 
•STATI 
XTCdD 
NE.C) 
96,204 
•MEAN 
XTOdD 
NE.C) 
96,206 
•MIN  P 
PERIO 


DEV,  IG 

OTPUT< 

DEV,  IT 

OUTPLT 

,ITEXT 

D«  ) 

EV,ITE 

LTPLT( 

,ITEXT 

1  ) 

EV,ITE 

UTPUT( 

,ITEXT 

STICS 

EV.ITE 

OUTPLT 

,ITEXT 

VALUE= 

EV,ITE 

OUTPLT 

,ITEXT 

ERIOC 

D=' ,F8 


CIR,5, IEP) 
10DIER,  "CGI1 
CIP,15,IEP) 


)     IEP,,D2, 


XT, 24, 3, 26, 1,2, IER) 

lOllIERt'PERICC* 

) 


XT, 24, 22, 75, 1,2, IER) 

lODIERt  «PTI^E« 

)NBLK 

ON    BLOCK  SI  2) 

XT, 24, 24,22, 1,2, IER) 

(101)IER,'TE1' 

)MFAN. SDEV 


2) 


)     M  N ,  M  A  X 

IN    ^ICSCSEC=I ,F8.3, 

.3) 
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20£ 


2  1C 


209 


42 
40 


41 
12C 


215 


216 


99E 


CAL 
IF( 

ek 

FCF 
(/L 
IF( 

enc 

FCF 
CAL 
IF( 
ENC 
FCF 
1MN 
CAL 
IF( 

cat 

IF( 
IF( 

CAL 
1F( 
CM 
CM 
CCN 
CEC 
FCF 
IF( 
CNA 
ENC 
FCP 
CAL 
IF( 
EK 
FCF 
CAL 
IF( 

ccn 


L  7E>TC 
1ER.NE. 
CDE(96, 
PAT(  «SE 
L  7EX7Q 
IER.NS. 
CDE(96, 
,VAT(  »TY 
L  TEXTO 
1ER.NE. 
CCE(96, 
NAT( »T0 
ANC  MA 
L  75>TO 
IER.NE. 
L  TEXTR 
IER.NE. 
KODCITD 
L  TEXT! 
PESP(l) 
N  =  MN 
>  =  NA> 
TIME 
CCE(  16, 
NAT(2F8 
CMIN.GE 
X2  =  C*-AX 
CCE(96, 
NA7(F8. 
L  7EXTC 
IER.NE. 
CCE(96, 
MT(F8. 
L  7E>70 
IER.NE. 

1INLE 


(IDEV»ITE 

0)  0U7PUT 
208, I7EX7 
7  MINI  S  M 
(ICEV, I7E 
0)  OUTPUT 
21CI7EX7 
PE  A  MIN, 
(IDEV,ITE 
C)CUTPUT( 
209,I7EX7 
LSE  S7A7 
X  7YPS  Jl 
(IDEV,I7E 
0)  OUTPUT 
<IDEV,777 
OOUTPUK 
IR(9) ,8). 
(IDEV,RES 
•NE.GH)  G 


120,RESF) 

•  0) 

•DMAXJ  GC 

-CM  IN 

215,I7EX7 

2> 

(ICEV,I7E 

C)0U7PU7( 

216,I7EX7 

2) 

(IDEV,I7E 

0)QUTPUTt 


X7 

(1 

) 

AX 

XT 

(1 

) 

A 

XT 

10 

) 

IS 

ST 

XT 

(1 

77 

10 

EQ 

o 


,24, 26, 24, 1,2, IER) 
01>IER,»TE2» 

VALUE  CF  Y  AXIS  IN 
,24,28,24, 1,2, IER) 
01  )ISR,,TE4» 


MCRSEC1  ) 


COMMA,  AND  NAX  IN  AN  F  FCRMATM 
,24, 29,24,1, 2, IER) 
1HER,  JPTYPE' 


TICAL  MIN 
A  C/R') 


ANC  VAX  AS  Y  AXIS 


A  L/K  '  ) 
,24, 22, 24, 1,2, IER) 
01)  IER,«T1« 
7778,1,20, 24,1, 2, IEP) 
DIER,  «TXTR« 
.ni  nn  m  z.n 


.0)  GO  TC  40 
4,0,9, IER) 
TC  41 


CMIN,DMAX 

TC  42 

)DMIN 

XT, 24, 20, 20, 1,2, IER) 
101JIER, 'LAS7* 
)CMAX 

XT, 24, 5, 20, 1,2,  IER) 
101JIER,  ♦LAST  +  I' 


ELFFERIN  4096  P0IN7S  CF  ARRAY  CCUN7  INTC  INCATA 

KEV=1 

C4LL    et'FF  ERIN  (1,1,  IND AT  A, NS ANP ,  INC) 
4     IF(  INC. EC. 1)    GO    TC    4 

IF(  (N8LK.LE.il) .CR. (NELK.GE. 14)  )    GC    TO    301 

CC    27    Ll=l,4081,16 

WPI7E(6i 106) ( INDA7A(L2),L2=L1,L1+15) 

FCFNAT(  IX, 16(16, 2X )  ) 

CCf^7INUE 

ccntinue 


1C6 

27 

8C1 


l«AKE    1NCA7A    A    REAL    ARRAY 

CC     10    I«1,NSAMP 
[A7A(I)=FLCA7(INDATA(I  ) )/4. 
1C    CCN7INUE 

IF( (KKK.EC.2).AN0.(NE^I.N£.LeLK) ) LR EC=N EW I-2*NSAM F 

ncf^;lize   (0  7C  n  tc  cisplay 

CC     12    I=1,LREC 

E  ATA  (I)  =  (DATA(I  )-CNIN}/CMAX2 
12    CCNTINUE 

FCFM    7f-E    AFRAY    INOATA    WITH    GRAPHICS    WCRDS 

L  =  l 
996    CCN7INLE 

CALL    CGINITIIDEV*  IGCIR,5,IER) 

1FUER.NE.0)    OUTPLT(lOl)     IER,'D1« 

7(1)=7INE 

CC     19    I  =1,iNDIS 

TINE=TINE+( ( (DATA U)*CMAX2)+DM IN)/ 1000000. ) 
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IF(I.EQ.(NCIS/2)  )T(2)=TIME 
IFU.EC.NDIS)    T(3)=TIMS 
19    CmiNlE 

ENCCCE(96,320,ITEXT)T(1),T(2)  ,T(3) 
3  20    FCFVAT<2  8XtF6,3,12X,F6.3,12X,F6.2) 

C/iLl    TE>T0<IDEV,ITEXT,24,22,l,l,i,IER) 
IF(IER.NE.C)0UTPUT(1C1)IER, •TTIVE1 
IhCATA(l)=IhEAD<0,101 
IhC4TM2l»XPACK(-«5t«5t0) 
INCATA(2)*IPACK<U./DIV)-.5,DATA(L),0) 
CC    *    I=2,NCIS 

INCATA(I*2)»IPACK((I/OIV)-«5,DATA(L+I-l)tll 
e    CCMINUE 

NVCFDsNCIS+2 


CF/SU    CN    SCREEN 


C4L 
1F( 

cn 

1F( 

N=h 

C/L 

IF( 

IF( 

IF( 

C-C 

L  =  L 

IF( 

KHK 

IF( 

GC 

S7C 

EUC 

>F, 


EC 


1€ 


'LCAC 


IER) 


•NS(  ICEVfISW,IERJ 

;.NE.G)_OUTPL7<101) 


L  FNS(ICEV,IS 
1ER.NE.C)  OUT 
LLS(  ISW,15).LT.C)  GC 


IER,»F1» 


GC  T1 


LLil  IbWt 15). L I  .0 ) 

LLS(  ISW, 16I.LT.01 

TO  80 

4NUIS 

L.lE.(NSAMP-j\IDIS  +  lJ  ) 

=  KKK-»1 

KKK.GT.2)NBLK=NeLK4l;GG 

70  998 

F 

PAP 


18 
995 


GO  TC  996 


TO  999 
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